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Abstract 

Applying Lie symmetry method, we find the most general Lie point 
symmetry groups of the radiation natural convection flow equation(RNC). 
Looking the adjoint representation of the obtained symmetry group on 
its Lie algebra, we will find the preliminary classification of its group- 
invariant solutions. 
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1 Introduction 

This paper can be viewed as a continuation of the work of paper [7] where 
the authors by using Lie group analysis for the PDE system corresponding 
to radiation effects on natural convection heat transfer past an inclined 
surface (RNC) could reduce this system to an ODE one and obtain nu- 
merical solutions by applying Runge-Kutta method. 

There are various techniques for finding solutions of differential equa- 
tions but most of them are useful for a few classes of equations and ap- 
plying these techniques for unknown equations is impossible. Fortunately 
symmetries of differential equations remove this problem and give exact 
solutions. 

Symmetry methods for differential equations were discussed by S.Lie 
at first. One of the most important Lie surveys was finding relationships 
between continuous transformation group and their infinitesimal genera- 
tors that allows us to reduce invariant conditions of differential equations 
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under group action, that arc complicated because of non linearity, to linear 
ones. 

The main results of applying Lie symmetry method for differential 
equations is finding group invariant solutions. A useful way for reducing 
equations is finding any subgroup of the symmetry group and writing 
invariant solutions with respect to this subalgebra. This reduced equation 
is of fewer variables and is easier to solve. In fact for many important 
equations arising from geometry and physics these invariant solutions are 
the only ones which can be studied thoroughly. 

This type of equations as (RNC) system plays an important role in 
engineering and industrial fields. In this paper Lie symmetry method is 
applied to find the most general Lie symmetry group and optimal sys- 
tem of (RNC). Finally, we gain group-invariant solutions of the reduced 
system. 



2 Symmetries of RNC System 

Consider the (RNC) system: 

du _|_ _ Q 
dx dy 
du du d^u . 

''d^^''d^ = W^ ^ ^ 

U—+V— - —(1 + 4R) — 
dx dy Pr dy"^ ' 

where Gr = ^^^^^ Taa)v Grashof number, Pr = ££ztL is the 

Prandtl number and R = — °° is the radiation parameter [7] . Now the 
infinitesimal method is implied as follow: 

The infinitesimal generator X on total space of the form: 

X = Y^C{x,u)d^i+Y,^'{x,u)d^, (2) 

1=1 i=i 

has the n*'' order prolongation Pr'^^X ([6], Th 4.16). Applying the fun- 
damental infinitesimal symmetry criterion ([6], Th 6.5) on X as follow: 

Pr("^X[A,(a;,M<"^)] =0, v = l,---,l, whenever A,(a;,M'"^) = 0. (3) 

obtains the infinitesimal generators of the symmetry group. 
The vector field associated with RNC system is of the form: 

X = ^i{x,y,u,v)dx + ^2{x,y,u,v)dy 

+(fii{x,y,u,v)du + (fi2{x,y,u,v)dv + ipz{x,y,u,v)dt- (4) 

Since in RNC system, second order derivatives are appeared, symmetry 
generators are obtained by applying (3) for the second prolongation of X. 
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The vector field X generates a one parameter symmetry group of RNC 
= (Ai, A2, A3, A4) if and only if equations (3) hold for v — 1,2,3,4. 
So symmetry group of RNC system is spanned by following infinitesimal 
generators: 

Xi = dx, X2 = dy, 

X3 = xdx+udu + tdt, X4 = 2xdx+ydy — vdv — 2tdt- 



The Lie table of Lie algebra Q for RNC system is given below: 





Xi 


X2 


X3 


X4 










Xi 


2X1 


X2 











X2 


Xz 













X4 


-2X1 


-X2 









Theorem 1. // gk{h) be the one parameter group generated by Xk, 
k = 1, - ■ ■ ,4, then 
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{x,y,u,v,t) 
{x,y,u,v,t) 

{x,y,u,v,t) 
{x,y,u,v,t) 



{h + x,y,u,v,t), 
{x,h + y,u,v,t), 

/ h h , h\ 

[xe ,y,ue ,iJ,£e ), 

/ '2h h —h I —2h\ 

(xe ,ye ,u,ve ,te ) 



(5) 



3 Optimal System 

The aim of Lie theory is classifying the invariant solutions and reducing 

equations. Some of Lie algebras contain different subalgebras, so classi- 
fying them plays an important role in transforming equations into easier 
ones. Therefore we are classifying these subalgebras up to adjoint repre- 
sentation and finding an optimal system of subalgebrgis instead of finding 
an optimal system of subgroups. 
The adjoint action is given by the Lie series 

Ad(exp(eXOXj) = Xj - e[Xi,Xj] + ^[X, [Xi,Xj]] - • • • , (6) 

where is the commutator for the Lie algebra, e is a parameter, 

and i,j = 1, •••,4 ([5], ch 3.3). 

Theorem 2. A one- dimensional optimal system of RNC results as fol- 
low: 

1) Xi, 4) X1+X2, 7) X3-X2, 

2) X2, 5) X2-X1, 8) X3+X4, 

3) Xa, 6) X2-I-X3, 9) X4-X3. 
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Proof: Consider the symmetry algebra Q of the RNC system with adjoint 
representation demonstrated in the below table: 





Xi 


X2 




X4 




Xi + sXs + 2eXi 


X2 


^3 


Xi 




Xi 


X2 + eX4 


X3 


Xi 


X3 




X2 


X-, 


Xi 






6-^X2 


X3 


Xi 



Let 

X = aiXi + 02X2 + 03X3 + aiXi. (7) 

In this stage our goal is to simplify as many of the coefficient Ui as possible. 
The following results are obtained from adjoint application on X: 

1. Let tti ^ 0. Scaling X if necessary, suppose that ai = 1. Then 

X' = 01X1 +02X2+ 03X3 +X4. (8) 

If we act on X' with Ad(exp(a3X3)), the coefficient of Xi can be 
vanished: 

X" = 02X2 +03X3+ X4. (9) 
By applying Ad(exp(a4X4)) on X", we will obtain X'" as follow: 

X"'=03X3+X4. (10) 

1-a. If 03 =^ 0, then the coefficient of X3 can be -1 or 1. Therefore 

every one-dimensional subalgebra generated by X in this case 
is equivalent to X3 + X4, X4 — X3. 

1- b. If 03 = 0, then every one-dimensional subalgebra generated by 

X in this case is equivalent to X4. 

2. The remaining one-dimensional subalgebras are spanned by vector 
fields of the form X with ai = 0. 

2- a. If as / 0, let 03 = 1. By acting of Ad(exp(o3X3)) on X, we 

will have 

X = a2X2+X3. (11) 

2-a-l. If 02 7^ 0, then the coefficient of X2 can be -1 or 1. So every 
one-dimensional subalgebra generated by X in this case is 
equivalent to X2 -I- X3 , X3 — X2 . 

2-a-2. If a2 = then every one-dimensional subalgebra generated 
by X3. 
2-b. Let ai = 0, 03 = 0. 

2-b-l. If 02 7^ 0, then we can make the coefficient of Xi either -1 
or 1 or 0. Every one-dimensional subalgebra generated by 
X is equivalent to Xi -f X2 , X2 — Xi , X2 . 

2-b-2. If 02 = 0, then every one-dimensional subalgebra generated 
by X is equivalent to Xi . □ 
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4 Invariant Solutions and Reduction 

Now for finding invariant solutions using characteristic method [2]: 

Qa \u=u(x,t)= X[u" ~ u°'(x,t)] \^^^(^^^t) a = l,---,M (12) 

where M is the number of dependant variables. Now invariant solutions 
for RNC result as follow: 

Consider Xi = dx- Applying characteristic method on X\ obtains: Qu = 

--^u{x,y) ,Qv = --§^v {x,y) ,Qt = --§^t{x,y). Solutions of this sys- 
tem are of the form: u = F\ {y) ,v = F\ (y) ,t = F\ {y). 

Consider X2 = dy. Applying characteristic method on X2 obtains: Qu = 
-^M [x, y),Qv = --§^v (x, y) ,Qt = -■§^t {x, y). Solutions of this system 
are of the form: w = Fi (a;) , v = Fi (a) ,t = F\ {x). 

Consider X3 = xdx + udu + tdt- Applying characteristic method on Xz 
obtains: Qu = u {x, y) - x-^u {x, y) ,Qv = -x-^v {x, y),Qt =t {x, y) - 
x-§^t{x,y). 

Solutions of this system are of the form: u = xFi (y) ,v = Fi (y) ,t = 
xFi (y). 

Consider X4 — 2xdx + ydy — vdv — 2tdt- Applying characteristic method 
on Xi obtains: Qu = -2x-§^u{x,y) - y-§j^u{x,y) ,Qv = -v{x,y) - 
2x-^v{x,y)-y-^v {x, y) ,Qt = -2t{x,y) - 2x-§^t {x, y) - y-§^t {x, y). 

Solutions of this system are of the form: u = Fi ^"^^ ,v = Fi (yyz) 

Consider X1+X2 = dx + dy. Applying characteristic method on Xi + X2 

obtains: Qu = -§^u{x,y) + -§^u{x,y) ,Qv = -§^v {x, y) + -§^v {x,y) ,Qt = 
■^t{x,y) + -^t{x,y). Solutions of this system are of the form: u = 
Fi {-x + y) , v = Fi {-x + y),t = Fi {-x + y). 

Consider X2 — Xi = dy — dx- Applying characteristic method on X2 — Xi 
obtains: Qu = -§^u {x, y) - ■§^u {x, y),Qu = -§^v {x, y) - ■§^v {x, y) ,Qt = 
-^t{x,y) — -^t{x,y). Solutions of this system are of the form: u = 

f'i {x + y) ,v = Fi {x + y) ,t ^ Fi {x + y). 

Consider X2 + X3 = xdx + dy + udu + tdt- Applying characteristic 
method on X2+X3 obtains: Qu = u {x, y) — x-^u {x, y) — -^u {x, y) ,Qv = 
-x-§^v {x, y) - {x, y),Qt = t {x, y) - x^t {x, y) - -^t {x, y). Solutions 
of this system are of the form: u = xFi (— In (x) + y) ,v = Fi (— In (x) + y) ,t = 
xFi{-ln{x) + y)- 

Consider X3 — X2 = xdx — dy + udu + tdt- Applying characteristic 
method on X3 — X2 obtains: Qu = u {x, y) — x-^u {x, y) + -^u {x, y) ,Qv = 
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-x-§^v {x,y) + -§:^v {x,y) ,Qt = t {x,y)- x-^t {x,y) + -^^t {x,y). Solutions 
of this system are of the form: u = xFi (In (x) +y) ,v = Fi (In (x) + y) ,t = 
xFi (In (a;) +y). 

Consider X3 + X4 — 3xdx + ydy + udu — vdv — tdt- Applying char- 
acteristic method on X3 + X4 obtains: Q„ = u{x,y) — 2>x-§^u{x,y) — 
y-^u {x, y) ,Qv = -v{x,y)-3x-^v (x, y) - y-§^v {x, y),Qt = -t {x, y) - 
^•'^'i^^i.^^y) ~ y'S^'t i^jy)- Solutions of this system are of the form: u = 
(#) ^,v = F^ (^) -i,^t = F, ^. 

Consider X4 — X3 = xdx + ydy — udu — vdv — Mdt- Applying charac- 
teristic method on X4 — X3 obtains: Q„ = —u(x,y) — x-g^u(x,y) — 
y-§^u{x,y) ,Qv = -V {x,y) - x-^v {x,y) - y-§^v ix,y) ,Qt = -3t{x,y) - 
x-^t{x,y) — y-^t{x,y). Solutions of this system are of the form: u = 
F,\l)x-\vJF,{l)x-\t = F,{r)x-\ 



5 Conclusion 



In this paper by applying infinitesimal symmetry methods, we find opti- 
mal system and finally could reduce the RNC system and find invariant 
solutions. 
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